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We perform the first study about the heavy-heavy-light quark potential in lattice QCD 
and a potential model. We find that the inter-two-quark confining force is reduced by valence 
quark motional efi'ects compared to the string tension. 

§1. Introduction 

There are a large number of studies about the interquark potential in hadrons. 
The interquark potential is mainly generated by gluon dynamics reflecting the SU(3) 
gauge symmetry. In particular, as many lattice QCD studies, the static and quenched 
approximation produces the only gluonic potential. However, realistic hadrons in- 
clude not only gluons but also quarks, and quarks affect the interquark potential. 
We can consider two types of the quark effects, dynamical quark effects and valence 
quark effects. An example of the dynamical quark effects is the string breaking, 
which is obtained by unquenched calculation. In this study, we focus on the mo- 
tional effects of finite-mass valence quarks, that is, valence quarks are not static but 
have a finite mass. 

To investigate such quark effects, we define the heavy-heavy-light {QQq) quark 
potential VQQg{R). It is defined as the energy of QQq systems in terms of the inter- 
heavy-quark distance R. The QQq system is constructed from two heavy quarks 
which are infinitely heavy and spatially fixed, and one light quark which is finite 
mass and moving around. Since this QQq potential effectively includes the light- 
quark motional effect, we can investigate the finite- mass valence quark effects to the 
interquark force. This QQq system is also related to the doubly charmed baryon, 
which has been recently discovered at the SELEX experiment. ^^'■^^ 

To calculate the QQq potential, we adopt the two approaches, SU(3) lattice 
QCD^) and a non-relativistic potential model. Lattice QCD is the first principle 
calculation based on QCD. The potential model, or the quark model, is one of the 
most successful analytical approaches for the low-energy hadron physics. 

§2. Lattice QCD approach 

In lattice QCD, the QQq potential is extracted from the expectation value of 
the QQq Wilson loop WQQg{R,T). As shown in Fig. [H the QQq Wilson loop is 
constructed from the heavy-quark trajectories and the light-quark propagator. The 
calculation is performed on 16^ quenched lattice with /3 = 6.0. For the statistical 
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Fig. 2. The Coulomb gauge fixed QQq potential 
Fig. 1. The QQq Wilson loop. The wavy line Vqqq in the lattice QCD approach. All the 
represents the light-quark propagator and scales are measured in the lattice unit a ~ 
the straight line the heavy-quark trajectory. 0.10 fm. 



error reduction, we apply the wall-to-wall propagator for the light-quark propagator, 
using the Coulomb gauge. Since the Coulomb gauge fixing can give a gauge variant 
contribution to the short-range physics, we investigate the long-range behavior of the 
potential, which is empirically unchanged by the Coulomb gauge fixing.^) Such a 
procedure surprisingly suppresses the statistical error. 

As the analogy of the QQ potential, the QQq potential is fitted with 

VQQqiR) = a,sR - ^ + Ceff. (2-1) 

As shown in Fig. [3] and Table HI this function is surprisingly suitable for the QQq 
potential. Compared to the string tension a of the static QQ or 3Q potential, the 
effective string tension a^fi is about 10-20% reduced, 

cjeff < 0- ~ 0.89 GeV/fm, (2-2) 

except for k = 0.1200 {Mq = 1.5 GeV). This result means that the inter-heavy-quark 
confining force is reduced due to the light-quark existence. 

As for the light-quark mass dependence, the effective string tension increases 
when the light-quark mass increases. In the infinite light-quark mass limit, the QQq 
system becomes the static 3Q system, and then a^s would approach a. 



§3. Potential model approach 



In the potential model, the QQq system is treated as a non-relativistic con- 
stituent quark system. The Hamiltonian of this system is written as 

where the subscripts 1, 2, and 3 mean the two heavy quarks and the light quark, 
respectively, and Mq is the light-quark constituent mass. For the three-quark inter- 
action, we adopt not the simple sum of the two-body force in ordinary quark models. 
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but the three-body force based on the static potential in quenched lattice QCD^^ 
as 

^(n, r2, r3) = aL^in - V , /'I , , (3-2) 

where a = 0.89 GeV/fm and A'^q = 0.13. The color flux-tube length Lmin is given as 
the length minimally connecting the three quarks. Once the heavy-quark coordinates 
are fixed, we can calculate the energy E{R) of the system only through the light- 
quark wave function. By minimizing E{R) in discretized space, we determine the 
ground-state light-quark wave function ipif^) and the ground-state QQq potential 
VQQq{R). For the energy variational calculation, we adopt the "renormalization 
group inspired variational calculation" .^^ 

An example of the light-quark spatial distribution is shown in Fig. [3l The QQq 
potential is fitted with Eq. ()2Tp in the range of i? < 1.2 fm, and the result is shown in 
Fig. [Hand Table HIl We can see the same behavior in lattice QCD, i.e., the reduction 
of the effective string tension and its Mg-dependence. To compare with the lattice 
QCD result, we calculate in the potential model with the same condition, Mg = 1 
GeV and R < 0.8 fm, and find a"eff — 0.76 GeV/fm. Thus, the calculation of the 
simple potential model almost reproduces the result of lattice QCD. 

In this potential model, we can calculate the expectation value (Lmin) of the 
flux-tube length from the obtained light-quark wave function. We find that (-Lmin) is 
almost a linear function of the inter-heavy-quark distance R, i.e., (Lmin) — ^'o + biR, 
and that the essential point is the relation between the flux-tube length and the 
inter-heavy-quark distance, as 

(crLmin) — {hia)R + boa ~ (JcsR + const. (3'3) 

Thus, (Tefi is reduced as cJcfr ^ bia, and bi means the reduction rate of aes- 

Note that the relation between the flux-tube length and the inter-heavy-quark 
distance can be generally more complicated, reflecting the light-quark dynamics. We 
confirm that, in the region of i? < 2.4 fm, the QQq potential deviates from Eq. (j2-ip 
and (Jeff is an increasing function of R. In the i? — > oo limit, R approaches (Lmin), 
so that cJeff would approach a. 



Table I. The lattice QCD result with the light-quark Table II. The potential model re- 
hopping parameter k of the 0(a)-improved Wilson suit. The error from the finite 
fermion.^' The constituent quark mass Mq is estimated volume and discretization is es- 
as the half of the p meson mass. timated within a few percent. 



K 


M, [GeV] 


a,ff [GcV/fm] 




x7^dof 


Mg [GeV] 


<7eff [GeV/fm] 




0.1380 


0.5 


0.73(3) 


0.13(1) 


1.16 


0.33 


0.73 


0.12 


0.1340 


0.7 


0.73(8) 


0.13(2) 


1.11 


0.50 


0.75 


0.12 


0.1300 


1.0 


0.75(8) 


0.13(2) 


1.18 


1.0 


0.80 


0.12 


0.1200 


1.5 


0.89(4) 


0.12(2) 


1.31 


2.0 


0.84 


0.12 
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§4. Conclusion 

From the two approaches, we can conclude that the effective string tension be- 
tween the two heavy quarks is reduced by the hght-quark effect, compared to the 
string tension. The reduction of the effective string tension originates from the geo- 
metrical relation between the inter-heavy-quark distance R and the flux-tube length 
Lmin- This relation is determined by the light-quark dynamics, and reflects the char- 
acteristic flux-tube structure in baryons/-''^^ Such a light-quark effect is expected 
to hold not only for QQq systems but also for ordinary baryons. In addition, also 
in multi-quark systems, the inter- two-quark confining force would be affected by the 
existence of other light quarks. 
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